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We have located the global minimum for all lead clusters with up to 160 atoms using a glue 
potential to model the interatomic interactions. The lowest-energy structures are not face-centred 
cubic as suggested previously. Rather, for TV < 40 the majority of structures are decahedral or 
hexagonal close-packed, and beyond this size the structures do not correspond to any of the structural 
forms commonly found in clusters. However, these latter clusters are not simply disordered. High 
symmetry, magic number clusters are still present, the most prominent of which is the 148-atom 
Did hexagonal barrel. We relate these structural preferences back to the form of the interactions. 

PACS numbers: 61.46.+w,36.40.Mr 
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I. INTRODUCTION 

The structure of a cluster is one of its primary proper- 
ties and one which has been intensely studied, experimen- 
tally and theoretically [1, 2]. However, there is still much 
to be learnt about the fundamentals of cluster structure 
and the possible structures that can be formed. For 
atomic clusters with pair interactions, it is relatively well- 
understood how the form of the potential determines the 
observed structure. For example, the effects of the width 
of the potential well [3, 4] and oscillations in the poten- 
tial [5-7] have been systematically studied. However, for 
the systems that are of most interest, the interatomic in- 
teractions are usually much more complex. In particular, 
metal clusters, which are of great technological relevance 
[8] , have a strong many-body character to their bonding. 

This presents a number of challenges to our under- 
standing of cluster structure in metals. First, there is the 
possibility that new types of structure could emerge as 
a result of many-body effects. Although, the structural 
types observed for pair potentials are also frequently ob- 
served for metals, e.g. the competition between icosahe- 
dral, decahedral and close-packed structures is also com- 
mon for metals [9], there are an increasing number of 
intriguing exceptions. One seemingly common feature 
for small metal clusters is to exhibit structures with no 
discernible order [10-14]. However, it might be that the 
disorder is a result of new structural principles that can- 
not be fully satisfied at the sizes considered (hence the 
disorder), but which could lead to novel high symmetry 
structures at certain magic sizes [15]. Indeed, there are 
general grounds to expect high symmetry structures to 
emerge irrespective of the potential [16]. Studies that 
just reoptimize known cluster structures will of course 
miss such new features, and so it is important that effi- 
cient global optimization algorithms are used. 
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Secondly, the many-body character makes it increas- 
ingly difficult to relate the observed structure back to the 
interactions, even when the assumed form for the many- 
body potential is relatively simple. There has been some 
interesting progress recently in this area, namely into the 
causes of the disordered structures [17], and the effect of 
the range of the attraction and repulsion on the compe- 
tition between icosahedral, decahedral and close-packed 
clusters [18] , but there is much still to be discovered. This 
task is particularly important because of the difficulty in 
producing good empirical metal potentials (it is not fea- 
sible to study the sizes in which we are interested in any 
other way). One needs this kind of physical insight to 
understand the strengths and deficiencies of a potential 
and how it could be improved. It would also help one to 
discriminate between different potentials that purport to 
model the same material but give rise to different struc- 
tures. 

Lead clusters illustrate some of these challenges. The 
first theoretical study on large clusters by Lim et al. us- 
ing a glue potential seemed to indicate that the most 
stable clusters at relatively small sizes (from at least 
N ~ 55) are face-centred cubic (fee) [19], the preferred 
bulk structure [20] . This conclusion was based on a com- 
parison of the energies of a series of Mackay icosahedra 
and fee cuboctahedra. It is quite unusual to see bulk 
structures already being favoured at such small sizes, but 
this finding was rationalized on the basis of the partic- 
ularly small value of 7, the ratio of the surface energies 
of the {100} and {111} faces [19]. Usually, the Mackay 
icosahedra have an energetic advantage at small sizes, 
because exclusively having {111} facets gives them an 
appreciably lower surface energy. Furthermore, these re- 
sults were not inconsistent with the experiments in the 
literature at that time, which were mass spectroscopic 
studies on very small lead clusters [21, 22], and electron 
diffraction experiments on very large clusters, which were 
identified as fee, but with possibly some vestiges of amor- 
phous structure [23]. 

This basic picture though has recently been challenged 
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by both experimental [24] and theoretical [25] results. 
Electron diffraction of clusters from 3 to 7 nm indicates 
that the largest clusters are dominated by decahedra, but 
for the smaller clusters it was not possible to obtain an 
adequate fit to the diffraction pattern, suggesting that 
alternative structural models need to be considered [24]. 

Simulations of the melting and freezing of large lead 
clusters (modelled by the glue potential used by Lim et al. 
[19]) unexpectedly revealed that for a certain size range 
(600 < N < 4000, at least) fee structures are not lowest 
in energy [25] . Instead, a new type of icosahedral struc- 
ture, which is more stable than the fee structures, spon- 
taneously formed both on freezing and on heating at tem- 
peratures just below that for melting. Similar structures 
had been previously seen in some simulations of large lead 
clusters but it was not recognised that they could be low- 
est in energy [26] . They resemble anti-Mackay icosahedra 
[27], which have a Mackay icosahedral core but with most 
of the outer layer in 'hexagonal close-packed' (hep) sur- 
face sites rather than the 'fee' sites that would continue 
the packing in the Mackay icosahedra. 

These results naturally raise intriguing questions about 
the structures of small lead clusters. What alternative 
structural models might be needed to understand the ex- 
perimental results? Do lead clusters, as modelled by the 
glue potential, really favour fee structures at small sizes? 
Lim et al. clearly showed that other standard forms were 
not lower in energy, so the structures would have to be 
somewhat unusual. Interestingly, a simple analysis us- 
ing macroscopic properties as inputs suggested that lead 
would be a particular likely candidate for disordered clus- 
ters to be low in energy [17]. 

Here, we address some of these issues by perform- 
ing global optimization for lead clusters with up to 160 
atoms. We pay particular attention to characterizing the 
structures of these clusters and to understanding why the 
potential favours the lowest-energy structures. 



II. METHODS 

To model the lead clusters we use a glue potential [28] 
of the form 



(i) 



where <j>{r) is a short-ranged pair potential, U{n) is a 
many-body glue function and m is a "generalized coor- 
dination number" for atom i. m is defined as 



n. 



(2) 



where p(r) is an "atomic density" function. These three 
functions have been fitted for lead using a variety of bulk 
and surface properties [19]. The use of surface energies is 
particularly important for the application of this poten- 
tial to model clusters. As well as clusters, this potential 




FIG. 1: The three functions that make up the lead glue po- 
tential: (a) (f>(r), p(r) and (b) U(n). The pair distances and 
rii values in the 13-atom decahedron are also plotted as im- 
pulses, with the heights proportional to the number that take 
that value. 



has been also used to model the surface reconstructions 
and pre- melting of low-index lead surfaces [29, 30], and 
lead nanowires [31, 32]. 

The choice of this empirical potential is motivated by 
the need for computational efficiency in order that global 
optimization is feasible for the sizes we consider here, and 
by our intention to compare with previous results. The 
use of ab initio electronic structure methods for lead is 
prohibitively expensive, especially as relativistic effects 
would need to be included to obtain reasonable results 
[20]. For example, sophisticated density functional calcu- 
lations are unable to reproduce the experimental surface 
energies and anisotropies [33, 34]. 

The functions U(n), <f)(r) and p(r) are displayed in Fig- 
ure 1. The pair potential has a very shallow well and so 
most of the binding energy comes from the glue term. 
The glue term has been chosen to have its minimum at 
n=12, consistent with the designation of n as an effec- 
tive coordination number. The form of p(r) is particu- 
larly significant. As p(r) decreases relatively slowly with 
increasing r beyond the minimum in the pair potential, 
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next-nearest neighbours make a significant contribution 
to n. Therefore, the difference in surface energies be- 
tween the {111} and {100} faces is small because, al- 
though an atom on a {111} face has fewer nearest neigh- 
bours, it has more next- nearest neigbours [19]. However, 
p(r) then decreases relatively rapidly to zero at the cutoff 
at r = 5.503A, which typically occurs between the second 
and third neighbour shells. 

For a pair potential the pair distances are the most im- 
portant quantities. For example, the lowest-energy struc- 
ture of a cluster involves a balance between maximization 
of the number of nearest neighbours, whilst minimizing 
the strain energy that results from nearest-neighbour pair 
distances deviating from the equilibrium pair value, r oq 
[3]. However, for a glue potential, such as the current 
one, where the main contribution to the energy is from 
the glue function, the most important quantities are the 
n,. Indeed, one of the key factors in generating a low- 
energy structure is to have the rii values as close as possi- 
ble to n oq , the value of n at the minimum of U. This can 
potentially lead to different ordering principles than for 
pair potentials. Only structures that have their nearest- 
neighbour pair distances close to r eq are generally com- 
petitive for pair potentials. However, this constraint is 
relaxed for glue potentials, and particularly when, as in 
the current case, p(r) initially falls off weakly with r. 

For the atoms on the surface of a cluster n, < n cq . 
Therefore, there will be a driving force for contraction 
of the surface to make the pair distances for the surface 
atoms smaller and hence their larger. At equilibrium 
the surface contraction will be balanced by the increase 
in energy due to the resulting compression of the cluster 
core. 

These considerations represent a particular problem for 
some of the usual forms for atomic clusters, such as the 
Mackay icosahedra and to a lesser extent decahedra. The 
inherent strain in these clusters results in pair distances 
between surface atoms that are longer than r eq , and so 
the compression needed to increase rn for the surface 
atoms is particularly large. Therefore, these traditional 
structural forms are expected to become increasingly dis- 
favoured by potentials for which the pair separation de- 
pends strongly on coordination number [17]. Instead, 
novel forms that are able to obtain large m values for 
the surface atoms, whilst not having too large an ener- 
getic penalty for compression of the cluster interior, could 
potentially be lowest in energy. 

The global optimization of the lead clusters was per- 
formed using the basin-hopping [35, 36] (or Monte Carlo 
minimization [37]) approach. This method has proved 
particularly successful in locating putative global minima 
for a wide variety of cluster systems [38] . The optimiza- 
tion task becomes rapidly more difficult with increasing 
N (e.g. the number of minima on the potential energy 
surface is thought to scale exponentially with N [39-41]) 
and so, of course, the possibility that we have not been 
able to obtain the true global minimum increases. How- 
ever, the structural principles and trends are clear from 
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FIG. 2: Energies of the putative global minima relative 
to Save, a four-parameter fit to these energies. Save = 
-2.0251iV + 1.6608iV 2/3 + 1.3662iV 1/3 - 0.8634 

our results. 



III. GLOBAL MINIMA FOR N < 160 

The energies and point groups for the putative global 
minima are given in Table I. Point files will be made 
available online at the Cambridge Cluster Database [38]. 
The energies of the global minima are represented in Fig- 
ure 2 in such a way that makes particularly stable clusters 
stand out. All clusters in the range 9 < N < 40 are de- 
picted in Figure 3 and a selection of particularly stable 
larger clusters in Figure 4. 

First, we will look at the global minima for N < 40 
in detail, before surveying more briefly the results for 
larger clusters. For N < 8 the clusters exhibit the same 
structures as typically seen for pair potentials. However, 
Pbg has a somewhat unusual form that can be described 
as two face-sharing octahedra, and so, as with Pbio, is 
the beginning of an hep cluster. Then for A?=ll and 12 
more open structures with three-fold axes of symmetry 
are preferred. 

Most of the global minima for 13 < N < 33 are dec- 
ahedral in origin. However, the growth sequence is not 
straightforward. The decahedra are generally asymmet- 
ric with the quasi-fivefold axis not passing through the 
centre of mass. So, although the growth sequence be- 
gins by adding atoms around the equator of the 13-atom 
Ino decahedron [42] , before this shell is completed, asym- 
metric decahedra with a longer quasi-fivefold axis be- 
come lower in energy, starting at N=21. Furthermore, 
sometimes part of the structure is distorted away from 
the ideal decahedral positions, e.g. at N=18, 19, 26 and 
32. There are also structures with two interpenetrating 
(Pb2o) and face-sharing (Pb 2 4) 13-atom decahedra, the 
latter with two additional shared capping atoms. It is 
noticeable that the decahedral global minima generally 
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TABLE I: Energies and point groups (PG) of the putative global minima. 



N 


PG 


Energy /eV 


N 


PG 


Energy /eV 


N 


PG 


Energy /eV 


N 


PG 


Energy /eV 


3 


D 3h 


-1.380851 


43 


C s 


-62.518056 


83 


Ci 


-131.453946 


123 


c 2 


-201.844885 


4 


T d 


-2.558548 


44 


C s 


-64.158589 


84 


D 2 


-133.335862 


124 


c 2 


-203.652142 


5 


D 3h 


-3.711742 


45 


C s 


-65.887013 


85 


Ci 


-134.977705 


125 


Ci 


-205.483134 


6 


o h 


-5.214277 


46 


Ci 


-67.545017 


86 


c 2 


-136.725287 


126 


Ci 


-207.149663 


7 


D 5h 


-6.342793 


47 


Ci 


-69.227454 


87 


Ci 


-138.426910 


127 


Ci 


-208.901679 


8 


Civ 


-7.665775 


48 


Ci 


-70.926964 


88 


Ci 


-140.139735 


128 


Ci 


-210.764087 


9 


D 3h 


-8.962242 


49 


c 2 


-72.624639 


89 


c 2t) 


-141.978803 


129 


Ci 


-212.620619 


10 


Civ 


-10.328111 


50 


Ci 


-74.303601 


90 


Cs 


-143.728686 


130 


Ci 


-214.378335 


11 


c 3v 


-11.771970 


51 


Ci 


-75.989247 


91 


C2v 


-145.469445 


131 


Ci 


-216.125616 


12 


D 3h 


-13.351511 


52 


Ci 


-77.748769 


92 


Ci 


-147.190877 


132 


Ci 


-217.911017 


13 


D 5h 


-15.060197 


53 


Ci 


-79.529727 


93 


Ci 


-149.017327 


133 


Ci 


-219.726303 


14 


Cl v 


-16.488673 


54 


Sio 


-81.438379 


94 


c 2t) 


-150.933479 


134 


Ci 


-221.493988 


15 


C2v 


-17.971698 


55 


Ci 


-83.050627 


95 


c 2t) 


-152.722173 


135 


Ci 


-223.228659 


16 


C2v 


-19.359118 


56 


Ci 


-84.670414 


96 


c. 


-154.419817 


136 


Ci 


-225.006412 


17 


c 3v 


-20.892141 


57 


Ci 


-86.343451 


97 


Ci 


-156.124859 


137 


Ci 


-226.856548 


18 


Ci 


-22.441282 


58 


Ci 


-88.160853 


98 


Ci 


-157.827739 


138 


c 2 


-228.680493 


19 


C2v 


-24.029140 


59 


Ci 


-89.961499 


99 


Ci 


-159.582628 


139 


Ci 


-230.562421 


20 


C2v 


-25.554526 


60 


Ci 


-91.725507 


100 


Ci 


-161.288978 


140 


Ci 


-232.461447 


21 


Civ 


-27.160557 


61 


Ci 


-93.369916 


101 


Ci 


-163.018692 


141 


Ci 


-234.347193 


22 


Ci 


-28.700367 


62 


Ci 


-95.068322 


102 


Ci 


-164.713475 


142 


Ci 


-236.226662 


23 


C2v 


-30.342369 


63 


Ci 


-96.835577 


103 


Ci 


-166.412388 


143 


Ci 


-238.007461 


24 


D 2h 


-31.834411 


64 


c 2 


-98.559680 


104 


Ci 


-168.225729 


144 


c 2 


-239.879381 


25 


C2v 


-33.394629 


65 


c 2 


-100.391525 


105 


Ci 


-170.097176 


145 


a 


-241.772812 


26 


C s 


-34.947504 


66 


s 4 


-102.045664 


106 


Ci 


-171.848423 


146 


c 2 


-243.665845 


27 


Civ 


-36.526823 


67 


Ci 


-103.803493 


107 


c. 


-173.738048 


147 


a 


-245.558654 


28 


Ci 


-38.036722 


68 


Ci 


-105.470725 


108 


c 2t) 


-175.653403 


148 




-247.451751 


29 


C s 


-39.653184 


69 


Ci 


-107.136603 


109 


C2v 


-177.373836 


149 


Ci 


-249.175286 


30 


C a 


-41.291166 


70 


Ci 


-108.921031 


110 


Cs 


-179.073251 


150 


a 


-250.900806 


31 


C2v 


-42.914946 


71 


c 2 


-110.659704 


111 


C2v 


-180.771206 


151 


Ci 


-252.632151 


32 


C2v 


-44.475555 


72 


Ci 


-112.360587 


112 


Cs 


-182.477331 


152 


c 2 


-254.384256 


33 


Ci 


-46.064990 


73 


Ci 


-114.048803 


113 


Cs 


-184.174089 


153 


Ci 


-256.115829 


34 


Ci 


-47.687973 


74 


Ci 


-115.771829 


114 


c 2 


-185.939956 


154 


c 2 


-257.846941 


35 


Cs 


-49.363201 


75 


Ci 


-117.552036 


115 


Ci 


-187.649103 


155 


Ci 


-259.573483 


36 


D 3d 


-51.115236 


76 


c 2 


-119.429639 


116 


Ci 


-189.351898 


156 


c 2 


-261.299215 


37 




-52.716480 


77 


Ci 


-121.236712 


117 


Ci 


-191.105099 


157 


Ci 


-263.019401 


38 


D 3d 


-54.314542 


78 


c 2 


-123.043586 


118 


Ci 


-192.855140 


158 


c 2 


-264.738883 


39 


Ci 


-55.857448 


79 


c 2 


-124.684422 


119 


Ci 


-194.607282 


159 


Ci 


-266.456122 


40 


c 2 


-57.474215 


80 


Ci 


-126.364903 


120 


Ci 


-196.361528 


160 


Ci 


-268.174733 


41 


Cs 


-59.130894 


81 


Ci 


-128.007123 


121 


Ci 


-198.211397 








42 


Cs 


-60.803249 


82 


1^3 


-129.829950 


122 


Ci 


-200.012938 









have a significant proportion of surface atoms in {100} 
type environments. This feature reflects the small energy 
difference between fee {111} and {100} faces noted ear- 
lier. For materials that more strongly favour {111} faces, 
the most stable decahedral form is usually a Marks dec- 
ahedron [43], because this structure maximizes the pro- 
portion of {111} faces, whilst retaining a relatively spher- 
ical shape. However, for lead the most stable decahedral 



clusters occur at iV=13 and 23 (Figure 2). 

The other set of ordered global minima found for 
N < 40 are the hep clusters at iV=35-38. Again, these 
structures are somewhat unexpected, particularly as the 
fee truncated octahedron is possible at 7V=38, but this 
is further evidence of a preference for structures with a 
significant proportion of {100}-like faces. 

Of the other global minima for < 40, Pbi7 is related 
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FIG. 3: The global minima for N < 40. Each cluster is labelled by the value of N. 



to the 11-atom global minimum, but it is hard to discern 
any overall order for those at N=28, 29 and 34. Pb3g and 
Pb4o are somewhat related to the preceding hep struc- 
tures, as is clear from the viewpoint chosen for Figure 3, 
but again there is little order apparent on the other side 
of the cluster. 

Beyond iV=38 none of the global minima that we have 
located can be assigned to any of the usual structural 
forms. However, it would be too simplistic just to char- 
acterize the clusters as disordered. From Table I we can 
see that high symmetry structures are still present. Fur- 
thermore, if the clusters were just disordered one would 
expect cluster properties to evolve fairly smoothly with 
size. However, it is clear from Figure 2 that there are 
"magic number" clusters that are particularly low in en- 
ergy. Unsurprisingly, these magic numbers often corre- 
spond to the high symmetry clusters. 

Although most of the clusters in this size range have no 
discernible overall order, there are common local surface 
motifs that are repeatedly visible. However, only at a rel- 
atively few sizes can these local preferences be assembled 
into a structure that has clear overall order. 

As for the smaller clusters the surface structures of 
the global minima reflect the particularly low value of 7. 
However, this does not lead to structures with large {100} 
faces, but rather to many surface atoms with {100}-likc 
environments. The surfaces are typically covered with a 
patchwork of squares and triangles. So on the flat regions 
of the surface it is common to see atoms surrounded by 
three triangles and two squares (there are two ways this 
can be achieved), rather than the six triangles or four 



squares, that are typical of {111} and {100} surfaces, 
respectively. 

Pb54 is somewhat related to the Mackay icosahedron. 
It has an uncentred 13-atom icosahedron at its centre 
and a clear five-fold axis of symmetry. Along this axis 
it looks similar to the the D§h structure that was found 
by Wolf and Landman as a low-energy isomer of the 55- 
atom Lennard- Jones clusters [44], and which is related 
to the icosahedron by a single rearrangement in which 
the structure is twisted around a five-fold axis. However, 
there is a canted arrangement of squares and triangles 
around the equator of the cluster. 

The axial configuration of Pbs4 seems to be quite a 
common motif, and similar patterns can be seen in one 
of the chosen views for ^=60, 78, 95 and 148, the last 
based on a six-fold rather than a five-fold symmetric ver- 
sion of the pattern, thus making the top surface flat, 
rather than pyramidal. As the size of these clusters in- 
creases the pattern is, of course, extended outwards. The 
resulting motif is clearest for the highly symmetric, 148- 
atom global minimum. 

Pbi48 is the most prominent magic number in this size 
range (Figure 2). In shape, the cluster is a hexagonal 
barrel. Although the outer surface has a clear sixfold 
symmetry, this is in fact broken by the octahedron at 
the centre of the cluster. 

As flagged in the introduction, an important aim of 
this paper is not only to characterize the global min- 
ima for this lead potential, but to understand how the 
observed structures relate back to the form of the po- 
tential. We start by examining the decahedral 13-atom 



FIG. 4: A selection of particularly stable global minima for N > 40. Each cluster is labelled by the value of N. For most of 
the clusters two perpendicular views of the structure are given. 



TABLE II: The contributions to the energy for a series of 
55-atom structures, namely the global minimum (Ci), the 
fee cuboctahedron (Oh), the Ino decahedron (Dsn) and the 
Mackay icosahedron (Ih)- The structures are denoted by their 
point group (PG). ) and (i?£ urf ) are the average atomic 

energies for atoms in the interior of the cluster and on the 
surface, respectively. All the energies are measured in eV. 



PG 


Energy 






E strain 




(Ui) 


(Ef» lk ) 




Cx 


-83.051 


-1.796 


216 


4.684 


-81.254 


9.072 


-1.885 


-1.405 


o h 


-82.559 


-4.338 


216 


2.142 


-78.220 


8.783 


-1.956 


-1.360 


D 5h 


-82.438 


-4.248 


219 


2.322 


-78.190 


8.800 


-1.933 


-1.365 


h 


-81.295 


-5.944 


234 


1.076 


-75.351 


8.500 


-1.946 


-1.333 



global minimum, for which the nj and rii values have 
been included in Figure 1 . It is noticeable that there is a 
significant dispersion of nearest-neighbour distances. In 
fact the longest distance is 11.7% longer than the short- 
est, which compares to a 2.3% difference for the same 
structure when optimized for the Lennard-Jones poten- 
tial. As expected from the discussion in Section II, the 
structure distorts to move the rn values as close to n cq as 
possible, rather than keeping the nearest-neighbour dis- 
tances near to the minimum of the pair potential. This 
is achieved by an expansion along the fivefold axis and a 
contraction of the equator of the cluster. This reduces the 
rii values for the two vertex atoms on the fivefold axis, but 
increases the rii values for the other ten surface atoms, 
while maintaining the n$ value for the central atom close 
to n oq . A similar anisotropy of the pair distances has 
previously been noted by Lim et at in their analysis of 
the lead cuboctahedra [19]; there is a greater contraction 



for the {100} faces of the cuboctahedra than the {111} 
faces because of the enhanced contribution to rii from 
next neighbours across the diagonals of the squares on 
the {100} faces. 

To understand why novel structural forms are observed 
for this lead potential, we take Pbss as an example and 
compare the contributions to the energy from a series 
of competing structures (Table II and Figure 5). The 
global minimum is based on the 54-atom structure illus- 
trated in Figure 4 but with an additional surface atom. 
Also possible at this size are a fee cuboctahedron, an Ino 
decahedron and a Mackay icosahedron. 

In Table II we have decomposed the pair energy into 
two components: 

-^pair — ^nn^ -^strain > (3) 

where n nn is the number of nearest neighbours, e is 
the depth of the pair potential, and -Entrain is the en- 
ergetic penalty for distances that deviate from r eq , the 
distance corresponding to the minimum of the pair po- 
tential. More formally, 

^strain = ^ e-<f>(nj), (4) 
i<j,rij<r 

where ro is a distance criterion that distinguishes nearest 
from next-nearest neighbours. For all the structures we 
consider, there is a clear separation between these coor- 
dination shells. 

As expected the pair energy only contributes a small 
fraction of the total energy. It is also noticeable that 
-^strain is of similar magnitude to -E pa ir ■ This is in marked 
contrast to what occurs for pair potentials, where mini- 
mization of the strain energy is a key element of a struc- 
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FIG. 5: A comparison of the properties of the 55-atom global 
minimum (gmin) to the fee cuboctahedron (Oh), the Ino dec- 
ahedron (Dsn) and the Mackay icosahedron (Ih), the same 
four clusters as in Table II. (a) The atomic energies, Ei, for 
each atom in the cluster. The atoms have been ranked by 
their distance from the centre of mass, with atom 1 being the 
closest to the centre, (b) (rii(rij < r)). (c) n pa i rs (r). 



ture's stability. Strain for the global minimum is partic- 
ularly large. 

Although the pair energy is small in magnitude, it is 
structure sensitive and so it can still determine the rela- 
tive stabilities of structures when the energies from the 
glue term, -Egiue, are similar. For example, the major 
component of the difference in energy between the 55- 
atom cuboctahedron and decahedron is the greater pair 
energy of the cuboctahedron. 

It is clear from Table II that the global minimum's 
stability is a result of its significantly lower glue energy, 
which is a result of the atoms being able to achieve n, 
values that are closer to the ideal value, n eq . However, 
this lower glue energy is partially offset by the higher pair 
energy resulting from the distortion of the pair distances 
that is necessary to achieve an increase in rij. 

If we look at the atomic contributions to the energy 
(Table II and Figure 5(a)) it is clear that the lower en- 
ergy results from a lower average energy for the surface 
atoms (particularly atoms 13-33 in Figure 5(a)), which 
outweighs the somewhat less favourable energies for the 
atoms in the interior of the cluster. 

It is also interesting to understand what structural fea- 
tures of the global minimum lead to the larger value of 
(rij). We analyse this in Figure 5(b) and (c), first by 
looking at the cumulative contribution to (m) from pairs 
with distances less than r: 



(ni(rij < r)) 



1 

N 



(5) 



It is particularly interesting to note that (nf(r)) for the 
global minimum only becomes largest beyond 4.926A. 
Therefore, although the contribution to (m) from dis- 
tances beyond this distance is small in magnitude, it is 
key in stabilizing the global minimum relative to the more 
conventional forms. 

We can analyse this further by considering n pa j rs (r), 
the number of pairs of atoms that are separated by less 
than r. It can be seen from Figure 5(c) that the num- 
ber of pair distances within the radius of the cutoff dis- 
tance for p is significantly larger for the global minimum 
than for the competing structures, which in turn corre- 
lates with the larger value of (n,). However, this is only 
true because the cutoff is located between the second and 
third coordination shell. The cuboctahedron, decahedron 
and icosahedron all have a relatively narrow distribution 
of nearest-neighbour distances, which leads to a clear dis- 
tinction between the second and third coordination shell. 
By contrast, the global minimum has a much more dis- 
perse nearest-neighbour shell and hence there is no clear 
distinction between a second and third neighbour shell. 
Instead, there is a steady increase in n pa i rs (r) beyond the 
start of the second neighbour shell. 

Although the above analysis has been presented for a 
single cluster, repeating this procedure for other sizes has 
confirmed the generality of the conclusions. 
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FIG. 6: A comparison of the energies of the global min- 
ima and the new icosahedra discovered by Hendy [25, 45], 
to series of high-symmetric structures, which include Mackay 
and anti-Mackay icosahedra, cuboctahedra, truncated octahe- 
dra (with regular hexagonal {111} faces), Ino decahedra and 
Marks decahedra. The energies are measured with respect to 
Scuboct = -2.0284JV + 1.7929iV 2/3 + 0.9714iV 1/3 - 0.6342. 



IV. CONCLUSIONS 

We have shown by the locating the global minima for 
small lead clusters interacting with a many-body poten- 
tial of the glue form that, contrary to the original conclu- 
sion of Lim et al., these clusters do not adopt fee geome- 
tries for N < 160. Instead, they form a series of novel 
structures that are a consequence of the many-body char- 
acter of the potential. These results naturally lead one to 
wonder at what size bulk-like fee structures will develop. 
To help us answer this question we have plotted in Fig- 
ure 6 the energies of the global minima, alongside those 
for a number of sequences of high-symmetry structures 
and the novel icosahedral forms that Hendy obtained by 
simulations of freezing [25] and by construction [45]. 

The figure confirms that the icosahedra and decahe- 
dra are always higher in energy than the best fee struc- 
tures (except at the smaller sizes considered in the last 
section) and that the fee truncated octahedra become 
slightly lower in energy than the cuboctahedra [19, 25]. 
More interestingly, we can clearly see that both the global 
minima we have found and the new icosahedra produced 
by Hendy [25, 45] are significantly lower in energy than 
the best fee clusters. Thus, extrapolations between these 
two sizes ranges [46] suggest that fee clusters might well 
not be the lowest in energy for intermediate sizes, and 
hence that fee clusters are not lowest in energy until at 
least N ~ 15 000 [45]. This is a surprising result, for 
although it is not uncommon to find small (N < 100) 
metal clusters that do not exhibit any of the usual clus- 
ter structures, it is unprecedented for this behaviour to 



persist up to such large sizes. 

The results are also relevant to the ongoing issue of dis- 
ordered metal clusters. Like recent theoretical results for 
gold [10-12], cadmium, zinc [14] and vanadium [13] many 
of our global minima do not fit with the fee, hep, decahe- 
dral and icosahedral structures that are often found for 
close-packed materials. However, to call these lead clus- 
ters disordered would be too strong because, although 
most of the clusters for N > 40 have no overall stuc- 
tural order, there are common local structural prefer- 
ences which at a few sizes result in highly symmetric 
ordered structures that are particularly stable. 

For pair potentials, clusters tend to retain a lattice 
structure away from the magic numbers, because it is un- 
favourable for the pair distances to deviate significantly 
from the equilibrium value. For example, most small 
Lennard-Jones clusters can be considered to be based 
upon Mackay icosahedra, either with an incomplete outer 
layer or covered by an ordered overlayer [47] . By contrast, 
for metal clusters the many-body character of the bond- 
ing can make it favourable to break the lattice structure 
away from the magic numbers so that the atoms can (in 
the language of the current potential) increase their ef- 
fective coordination numbers (rij's). A structure with no 
overall order results. Therefore, if one only examines a 
few cluster sizes the presence of particularly stable or- 
dered structures may be missed. 

From our analysis of the energetics of the competing 
structural forms we have seen that the shape of p(r), 
in particular the shoulder and the position of the cut- 
off, is key to the stability of the novel structures that 
we find to be lowest in energy. This dependence on the 
cutoff is somewhat worrying both because it is a rather 
long-range feature of the potential and because its po- 
sition is not physically motivated, but chosen more for 
computational convenience. These results illustrate how 
sensitively cluster structure depends on the potential; the 
correct determination of the relative energies of compet- 
ing clusters is a stringent test of any potential [48] . 

Our results also help us to understand the structures 
exhibited by lead nanowires [31] modelled using the same 
potential. Gulscrcn et al. were surprised at the appar- 
ent contradiction between the non-fee character of their 
nanowires and the cluster results of Lim et al. [19] , and so 
suggested a number of reasons for the differing structural 
tendencies. However, our results show that the non-fee 
character is common to both systems. 
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